Anomalous Diffusion of Inertial, Weakly Damped Particles 
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The anomalous (i.e. non-Gaussian) dynamics of particles subject to a deterministic acceleration 
and a series of 'random kicks' is studied. Based on an extension of the concept of continuous 
time random walks to position-velocity space, a new fractional equation of the Kramers-Fokker- 
Planck type is derived. The associated collision operator necessarily involves a fractional substantial 
derivative, representing important nonlocal couplings in time and space. For the force-free case, a 
closed solution is found and discussed. 
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For more than 100 years, the quest to understand 
the dynamics of tracer particles in random environments 
has both challenged and inspired countless scientists. 
Many fundamental questions about the dynamical as- 
pects of statistical physics have been raised in this con- 
text, and most of them are still unanswered today. Given 
the fact that this nonequilibrium phenomenon plays a 
prominent role in such diverse areas as condensed mat- 
ter physics, biophysics, physical chemistry, astrophysics, 
plasma physics, and turbulence research, it is not at all 
surprising that its theoretical description has continued 
to attract a lot of attention over all these years. 

In the more recent past, a strong emphasis in this area 
of research has been on the investigation of the origin 
and the consequences of anomalous, i.e. non-Gaussian, 
diffusion (see, e.g., [ll|j,|3j). While most publications on 
this topic are based on a random- walk type description in 
real space, this Letter employs a phase (position-velocity) 
space approach, rigorously deriving and analyzing a new 
fractional equation of the Kramers-Fokker-Planck (KFP) 
type. For later reference, we mention here that the usual 
KFP equation is given by [J, ^ 



d_ 
di 



u ■ V^ + A(x) • V„ 



/(x,u,t) = /:Fp/(x,u,t) 



(1) 

where we have introduced the Fokker-Planck collision op- 
erator 

/:fp/ = 7V„-(u/) + kA„/. (2) 

As Kramers 6] has shown, Eq. QJ determines the time 
evolution of the joint position- velocity probability distri- 
bution /(x, u, t) of a damped particle with inertia in a 
force field under the influence of an additional fluctuating 
force with the usual white noise statistics: 



i(t) = u(t) , u(t) = A(x) - 7u(i) + r(t) (3) 



with 



In 1965, MontroU and Weiss '7!| introduced a different 
class of stochastic processes denoted as continuous time 
random walks (CTRWs). - For a recent, comprehensive 
review, the reader is referred to Ref. 3]. - Conventional 
CTRWs are linked to Eq. (Q in the limit of large damping 
7 and vanishing acceleration A(x). Here, the dynamics 
is described by the Langevin equation x(t) = T{t)/^. In 
practice, however, many systems in physics, chemistry, 
and biology will not satisfy this condition. Thus we are 
led to consider the CTRW analog of Eq. (j3l . 

The dynamical model we propose contains both deter- 
ministic and stochastic elements. The motion of an in- 
dividual particle is assumed to be governed by the time 
evolution equations 



x(i) = u{t) , u(i) = A(x) . 



(5) 



However, from time to time, the particle is subject to a 
'random kick' which changes its velocity abruptly. This 
kind of dynamics can be described by means of a CTRW- 
like approach in phase space. Let us consider a particle 
which at time t' is located in the volume element dx' 
about x', and in the time interval [t' , t' + dt'] changes its 
velocity to a new value which lies in the velocity space 
element du' about u'. The probablity for such a process 
shall be denoted by 



r]{x',u',t')dx'du'dt' . 



(6) 



(r^{t)Tj{t')) = 2KS^jSit-t') 



(4) 



After a (random) time period t — t — t', this particle will 
undergo a further transition to a state with the veloc- 
ity u at the position x. The corresponding conditional 
probability shall be called 

^(x, u, t;x', u') dxdudr . (7) 

We assume that this quantity can be written in the form 

C(x, u, r; x', u') = W{t) S{x - X^(x', u')) 

X / du" F(u; u") Siu" - U,(x', u')) • (8) 



Here, W{t) dr describes the probability that a transition 
occurs in the time interval [t^t + dr], and F(u;u")(iu 
gives the probability that the particle's velocity will end 
up in the velocity space element du about u. Moreover, 



x = X^(x',u'), u = U^(x',u') 



(9) 



is the solution (at time t) of Eq. |SJ with initial values 
(at time t') x' and u'. Assuming that the distribution 
function ^(x, u, r; x', u') is statistically independent from 
the particle's path and using Eq. JHl, we can relate the 
quantities 77(x, u, t) and ??(x', u', t') via the equation 

77(x,u,t) - /o(x,u)^(t) = 

= j dt' f du' /"dx'^(x,u,t-i';x',u')?7(x',u',t') = 

= f dt'W{t-t') f dWF{u;u')P''''T]{K,W,t'). (10) 

Here, /o(x, u) characterizes the initial condition, and 
the deterministic evolution is described by the Perron- 
Frobenius operator 

7'*^*'r,(x,u,t') = e-(*"*'^["-^^+^(")-^"l?7(x,u,t') = 

= f dW f dx'6{x-Xrix',u'))6{u-Vrix',u')) 
X7,(x',u',t'). (11) 

Note that, alternatively, one can write 

P*-*',(x, u, = ^(X.-nx,u)U-|(x,u),tO ^^2^ 

\Jt\X, ujl 

where (X^^,U^^) denotes the inverse of the mapping 
© , and Jr is the corresponding Jacobian. [For a Hamil- 
tonian dynamical system, Jt equals unity.] 

Having established an integral equation which deter- 
mines the time evolution of ?7(x, u,t), we are now inter- 
ested in the joint position-velocity distribution function 
/(x, u, t) which is defined as 

/(x,u,i)= / dt'w{t-t')r*'*'ri{x,u,t'). (13) 

Here, u;(t) denotes the probability that no random kick 
occurs within the time interval r. For this quantity we 
have the obvious relationship 



w{t) = 1-1 dtW{t) 



(14) 



In order to derive an equation describing the time evolu- 
tion of /(XjUjt), we first write down the Laplace trans- 
forms of Eqs. l(TU|l and ((T^ . 

77(x, u, A) = /o(x, u) + y du' F(u; u') 

X W{X + u' ■y^+ A(x) • V„/) ?7(x, u', A) (15) 



and 

/(x, u. A) = u;(A + u • V, + A(x) • VJ ?7(x, u, A) . (16) 
Since in Laplace space, Eq. ![T^ reads 

1 - W{X) 



w{X) 
Eq. Hlt)|) can be rewritten as 



A 



(17) 



iX + u-y^ + A(x) • V:,) /(x, u. A) = 

= [1 - W^(A + u • V, + A(x) • V„)] ^(x, u. A) . (18) 

Eqs. (HSJ and (HHJ then yield 

(A + u • V:, + A(x) • V„) /(x, u. A) = 

= /o(x, u) + /" du' F(u; u') $(A + u' ■ V, + A(x) • V„0 

/(x, u', X)~W{X + u-y^+ A(x) • V„) r,(x, u. A) (19) 

where we have introduced the quantity 

XW{X) l-Xw{X) 



$(A) 



wiX) 



1 - WiX) 

Expressing /o in terms of / as 

/o(x,u) =/(x,u,t = 0) 

and using the identity 

WiX -I- u • V^ + A(x) • V„) ?7(x, u, A) = 
= $(A + u . V, + A(x) • V„) /(x, u. A) 



(20) 



(21) 



(22) 



which follows from Eqs. H18() and (|20|l . the Laplace inver- 
sion of Eq. ifT^ yields 



d 



— +u-\7^ + A(x) • V 



/(x,u,t)- 

/ dt'$(i-i') f du' Fiu;u')V'''' fix,u',t') - 
[ dt'$(t-i')^*'*7(x,u,i')- (23) 



Thus the time evolution of /(x, u, t) is given by this mas- 
ter equation. In the language of kinetic theory, Eq. l|^ 
can be interpreted as follows. The left-hand side de- 
scribes a system of particles subject to Eq. ©. The 
right-hand side represents a collision operator which con- 
sists of a source and a sink. The phase space density of 
particles at (x, u) is increased at time t by particles start- 
ing from X~-'^(x, u) at time t' with a velocity U~^(x, u) 
and making a transition to the velocity u at time t and 
position X. /(x, u,t) is decreased, on the other hand, by 
particles making a transition from the velocity u to some 



other velocity. In this context, it should be pointed out 
that 

f dt'^t^t') /"du"F(u";u)P*'*'/(x,u,i') = 

= / dt'$(i-i')7'*'*'/(x,u,t') (24) 

Jo 



due to the constraint 



fdu"F{u";u)^l, (25) 

and that the Laplace inversion of Eq. (|20|l yields 
dw{t) 



dt 



= - / dt'<^{t-t')w{t') 



(26) 



which identifies <i>(i) as a kernel determining w{t). Also 
note that 7-"*'* /(x, u, t') is simply a solution of the coUi- 
sionless equation 



d 

— + U • V:, + A(x) • V„ 



/(x,u,i)=0. 



(27) 



Obviously, the above collision operator is highly nonlo- 
cal in space and time, in stark contrast to virtually all 
expressions commonly used in the kinetic theory of gases 
or plasmas. This nonlocality can also be viewed as a 
retardation effect and is closely linked to the fact that 
Eq. (|23|l is invariant with respect to Galilean transfor- 
mations, provided f (u; u') depends only on the velocity 
difference u — u'. [Galilean invariance implies that to- 
gether with /(x, u, t), every /(x — ct, u -I- c, t) with c — 
const is a solution of Eq. (|23|l .] 

For concreteness, we now consider the case where 
F(u; u') can be represented by the Gaussian 



(u;u').r-^^ 



3/2 



^^"'"') = ^4;^; ^"p 



(u - u' + -iW/kf 
4k/A 



(28) 



which satisfies Eq. H25|l . In the limit A ^ oo, one obtains 

j du' F(u; u') 5(u') - 5(u) = A-'Cppgiu) (29) 

to leading order in A~^. Removing the prefactor A~^ by 
means of the parameter redefinition 

7/A^7, k/A^k, (30) 

Eq. H23f) then takes the form 

' d 

— +u-V,+A(x)-V„ /(x,u,t) = 
ot J 

= £fp / dt' $(t - t') P*'*7(x, u, t') . (31) 

This equation, a rigorous result based on a consideration 
of continuous time random walks in phase space, is a 



nontrivial extension of the usual KFP equation, Eq. 1^, 
which is recovered for ^(t) = (5(r) [corresponding to 
w{t) — exp(— r)]. While generalized KFP equations with 
memory effects have been considered before by many au- 
thors, the retardation effect has never been included. The 
latter is important, however, in order to ensure Galilean 
invariance, as pointed out before. Moreover, retardation 
enters quite naturally in the present CTRW framework 
and its physical origin is evident, given the mixed nature 
of the underlying physical process - a particle being sub- 
ject to a deterministic acceleration and a series of random 
kicks. 

Let us now consider the important case that the mem- 
ory kernel <I> in Eq. (|31|l exhibits a power law tail, i.e. 



$(i - t') 



1 d 



1 



1 



r{S) dt {t - t'y- 



(t-f) 



i\2-5 



(32) 



in the long-time limit, where the prefactor is given in 
terms of the gamma function F. A regularization proce- 
dure J8!| then leads to the correspondence 



dt' $(< - t') H{t') -^ Dl-^H{t) 



(33) 



.1-5 



where the so-called Riemann-Liouville operator D^ is 
defined via 



Dr'H{t) 



1 d 



/■* dt' 



F(<5) dt 

for < ^ < 1. Consequently, Eq. H2t)|l turns into 
dw{t) 



dt 



+ D^-''w{t) = 0. 



(35) 



In the limit (5^1, one gets I?^^ w(t) — > w{t), and 
Eq. (^5)1 has the solution w{t) = exp(— i), describing 
classical relaxation behavior. However, for S < 1 one ob- 
tains fractional equations instead. A careful investigation 
of the regularization process actually shows that, in the 
present case, it is appropriate to introduce a novel type of 
fractional derivative - which might be called "fractional 
substantial derivative" - via 



i?r^ff(x,u,i) 



dt' 



1 



ns) 



d_ 
di 



u-Vx+ A(x) • Vt, 



it-t') 



Al-<5 



-(*-*''["-^^-+^W-^"li/(x,u,t'). (36) 



We note that, alternatively, one may define this fractional 
substantial derivative by means of its representation in 
Laplace space, 



V 



1-5 



^ [A + u ■ Vj; + A(x) ■ Vi 



n-s 



(37) 



Using this definition, one arrives at 

' d 

— +u-Wa,+ A(x) ■ V„ 

= CFpVl-'f{^,u,t), 



/(x,u,t) 



(38) 



a fractional generalization of the standard KFP equation, 
Eq. 1^. The latter corresponds to the special case 6 = 1. 
Moreover, it might be worth pointing out that the oper- 
ators £fp and 'DI~ do not commute. Eq. (|38|l differs 
significantly from the fractional KFP equation proposed 
by Barkai and Silbey [9j in that it includes retardation ef- 
fects through the fractional substantial derivative V}' . 
These nonlocal couplings in time and space are crucial, 
as will be shown next. 

For simplicity, let us focus on the case A(x) = in 
which the position-velocity distribution with the initial 
condition /(x, u, f = 0) = (5(x) (5(u) can be written as a 
superposition of Gaussians with varying variances: 



and normal diffusion in the long-time limit, while for 
(5 < 1, the corresponding probability distributions are 
non-Gaussian. In the latter case, the terms proprotional 
to (1 — (5) will dominate for very long times, leading to 
ballistic motion - independent of the values of 7 and k 
as long as both are non-zero. These terms have their 
origin in the collision term's inherent memory. Conse- 
quently, retardation effects are observed to alter the dif- 
fusion properties of the system even qualitatively. 



/(x,u,t) = da dp d'yw{a,P,j,t)[{2nf det Ay 



1/2 



X exp 



-— u ~qi2U- (x-ut) 



922 



(x-ut)2 . (39) 



Here, A denotes the matrix with the elements An = a, 
^22 = 7, and A12 = A21 = (3; moreover, qij are the 
elements of the inverse matrix Q = A~^. The function 
w{a, P, 7, t) can be viewed as a time-dependent variance 
distribution. It obeys a partial differential equation with 
temporal memory which can be solved analytically. For 
details, the reader is referred to Ref. [lOj . Eq. (|39|l can be 
used to calculate the long-time behavior of the second- 
order moments of the distribution function /(x, u, i). 
Here, it proves useful to distinguish two distinct situa- 
tions, 7 = and 7 7^ 0. For 7 = 0, we obtain 



lim {u{ty 

t— 5-00 



lim {x{t)u{t)) 



lim {x{tf) 

t — *-oo 



2k 



Til + S) 
2k 


r{2 + s) 

4k 



t^ 



fl+5 



i2+(5 



r(3 + 5) 



(40) 
(41) 
(42) 



In summary, using a CTRW-like approach, we were 
able to rigorously derive a new fractional equation of the 
Kramers-Fokker-Planck type for particles subject to a 
' deterministic acceleration and a series of random kicks. 
In this context, we found that the generalized Fokker- 
Planck collision operator necessarily involves a fractional 
substantial derivative, representing important nonlocal 
couplings in time and space. For the force-free case, a 
closed solution was given in terms of a superposition of 
Gaussian distributions with varying variances. This re- 
sult was then used to calculate the long-time behavior of 
the second-order moments, revealing the ballistic nature 
of the associated diffusion processes for any < 6 < 1 
as long as 7 and k are both non-zero. This universality 
could be identified as a direct consequence of retardation 
effects which are expressed mathematically by means of 
the fractional substantial derivative. In conclusion, the 
present work may be considered as a useful starting point 
for further investigations of anomalous diffusion of iner- 
tial, weakly damped particles. 



in the one-dimensional case.jlOJ The extension to multi- 
ple dimensions is straightforward and only modifies the 
prefactors in the above equations. For 5 > 0, the system 
exhibits superballistic diffusion; for (5 = 1, it reduces to 
the so-called Obukhov model. [IJ We would like to point 
out that for 7 = 0, the second-order moments are not 
affected by retardation effects, unlike the higher-order 
moments. For 7 7^ 0, on the other hand, one gets [l^ 



lim (uitf) = - 



(43) 



lim (x(i)u(i)) = j..^ ^, ^t'-' + {l-5)-tlU) 
t^oo 1 (2 — 0) 7^ 7 



lim {x{t)^) 



r(3 - 6) 7 



^ t'-' + {1-S)- t\A5) 



K 

7 



For 6 — 1, Eq. (|38(l turns into the usual KFP equation, 
Eq. (^, which is known to exhibit Gaussian statistics 
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